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, I. INTRODUCTION AND OVERVIEW 

o 



o 



It can be argued that, aside from theories with infinite number of particle types such as string theory, quantum 
field theory is the only way to reconcile the principles of quantum mechanics with those of special relativity [|. In 
, quantum field theory, vanishing correlations for space-like separated operators are ensured^]; whereas most efforts 
in quantum computation have so far relied upon non-relativistic quantum mechanics which is not fully compati- 
ble with Lorentz invariance and the causal structure of space-time. Recently however, several groups(see, for in- 
stance, Refs .|3.l3H3.l5Llr3.F^.lcll^. ITfl ITU ll2L IT^I ) have focused their investigations on relativistic effects in quantum 
information science. The issues include Lorentz invariance of entanglement, the behaviour of Einstein-Podolsky- 
Rosen-Bell states in different incrtial frames, and possible modifications to the degree of Bell inequality violations 
for moving observers^]. These relativistic effects may alter the efficiency of eavesdropper detection in quantum 
cryptography^^ and compromise the security of quantum protocols. It is also expected that future applications 
in quantum teleportationflaj . entanglement-enhanced communication [bj. high-precision quantum clock synchroniza- 
tion based on shared entanglement, and quantum-enhanced positioning [18j will also require relativistic treatments of 
quantum systems; and in particular, the careful analysis of the properties of entangled particles under Lorentz trans- 
formations and the construction of meaningful measures of entanglement. In this article, we study the behaviour of 
Bell states under Lorentz transformations and consider von Neumann entropy as a Lorentz invariant characterization 
of entanglement. 

We compute, for massive particles, the explicit Wigner rotations of one-particle states; and the explicit generators 
of the unitary representation of the Lorentz grou p. U nitary representations of the Poincare group acting on physical 
states are founded upon Wigner's seminal work[l!|; but the explicit expressions of Wigner rotations for massive 
particles have been computed, with some difficulty by direct matrix multiplication, only for rotations and boosts 
considered separately!^. Our derivation, carried out in Section II, is somewhat simpler, and it permits the explicit 
general result for arbitrary infinitesimal Lorentz transformations to be stated as in Ea. l|2.7|l . Moreover, with the 
infinitesimal Wigner angle at hand, the explicit infinite-dimensional Hermitian generators of the unitary representation 
can be worked out, as in Section III. Wigner rotation for a finite general Lorentz transformation is slightly more 
complicated, but the explicit form is also listed in Appendix A. 

As basic entangled states, Bell states figure prominently in the literature on quantum information science and 
in quantum computational schemes. Their behaviour under Lorentz transformations is therefore of interest and 
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importance. In this article we focus on Bell states of two spin | massive fermions. Under Lorentz transformations, 
each one-particle state of the entangled pair undergoes an SU(2) Wigner rotation. Since [SU{2) x 5J7(2)]/Z 2 is 
isomorphic to 50(4), which contains 3-dimensional rotations SO(3) as a sub-group; in group theoretical terms the 
entangled system transforms as 2(g) 2 = 4 = 103, in which the final step denotes the behaviour under SO(3). Thus by 
forming a 4-vector of SO(A) out of the rotational singlet and triplet Bell states, the complete, and explicit behaviour 
given the Wigner angle, of these Bell states under Lorentz transformations can be succinctly stated, as in Ea. H4.4JI . 
We therefore advocate for Bell states the convention in Ea. H4.2J) . Details of the setup are presented in Section IV. In 
Section V, we introduce reduced density matrices, defined through Yang's prescription |2l| . for systems of identical 
particles; and analyze their behaviour under Lorentz transformations. Without taking Lorentz symm etry into account, 
quantum correlations and entanglement of identical fermions have also been studied in Refs. |2^l23Ll24l l25|. Here we 
are able to show that the von Neumann entropy of each of the constructed reduced density matrix is Lorentz invariant; 
and we illustrate the usefulness of the von Neumann entropy as a Lorentz-invariant measure of entanglement with 
a worked example comparing unentangled and Bell states. In the ensuing subsection, we present a relation between 
generalized zeta functions and von Neumann entropy which may be useful for the regularization of the entropy of 
infinite-dimensional density matrices. Further comments and conclusions are presented in the final section. 

As far as it is convenient to do so, we shall follow the conventions and normalizations in Weinberg's tome.l], with 
if™ = diag(—l, +1, +1, +1). Space-time Lorentz indices are denoted by Greek letters, spatial indices by Latin letters; 
and summation over repeated discrete index is assumed. Our computations shall concentrate, for convenience, only on 
the homogeneous Lorentz group since translations can be incorporated rather readilyjlj once the behaviour of fields 
under the homogeneous group has been worked out. 



II. WIGNER ROTATIONS 



In quantum field theory, one-particle states \p, s) are classified by eigenvalues of the Casimir invariants of the 
Poincare group. Any value of the momentum, p 11 can be reached by Lorentz transformation L(p) on a standard fc M 
for which p^ — L^(p)k u \ and the states can be defined as 

\p,s)=J^U(L(p))\k,s), (2.1) 
V P 

with P M |p, s) = p^\p, s). It follows(see, for instance, Ref.0) that the effect of an arbitrary Lorentz transformation A 
unitarily implemented as U(A) on one-particle states is 



\p, s)' = U(A)\p, s) = J Q$-D*,(W(A,p))\Ap, s'), (2.2) 

and 

W(A,p) = L-\Ap)AL{p) (2.3) 

is a Wigner transformation which leaves A; M invariant |19j, and D(W) represents its action on the state (summation 
convention over the repeated index s' is assumed). The explicit form of L(p) is dependent on the class of the four- 
momenta. For massive particles, p^p^ — —m 2 < 0, and a convenient choice for the standard vector is k^ = (m, 0). 
It is then obvious that the set of Wigner transformations leaving k^ unchanged is just the rotation group SO(2>). 
Furthermore, L(p) can then be taken as the pure Lorentz boost 

L\{p) = cosh X (2.4) 
L\(p) = L i (p)=p iS inh X (2-5) 
L l } (p) = # j + (coshx -l)p%- (2.6) 

with tanhx = , ^ In this parametrization, L(p) — exp(— «xp • K) where K l — A/° l is the boost generator 

vipi 2+m2 

We may rely upon the analytic nature of Lie groups for the computation of the Wigner angle, and the corresponding 
infinitesimal Wigner rotation can be evaluated by Taylor expansion. Details of the rest of the computations are 
delegated to Appendix A. The end result is that the infinitesimal Wigner rotation of a massive particle is 



W(A,p) =I+| 



M ij (2.7) 



1 

p" 

1 + iOw-J, (2. 
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with the infinitesimal Wigner angle denoted as 

e w = e-^-^ = e + <f> 1 . (2.9) 

p + 111 

This agrees with the results of Ref. [2(j when rotations and boosts are considered separately. It should be noted that the 
generators of the Wigner rotations are indeed J; = \tijkM^ k , but the complete Wigner angle receives contributions 
from both the boost and rotation parameters, t 1 = uj 1 o and 6i = ^eijkUJ^ respectively, of A(u). In the absence 
of boosts, Wigner rotations are degenerate with ordinary rotations i.e. Bw = 0. Although arbitrary finite Lorcntz 
transformations of A(w) = exp(|w Q( aM a/3 ) can be evaluated as limjv->oo[I+f < ^-M a ^] N , it is not simple to express the 
corresponding finite Wigner rotations in closed form via products of infinitesimal rotations, essentially because Wigner 
rotations are also functions of the momenta. We may however consider a general Lorentz transformation relating two 
frames as a product of a pure boost in an arbitrary direction, L(a) = exp(— iot ■ K), followed by an arbitrary rotation 
R(ip). Using the multiplication rule for Wigner transformations, fEa. (|A7|) in Appendix A), the closed form Wigner 
rotation can be expressed as in Eqs.(A16)-(A17). Note also that for the special case of U(A) = U(L(p)) acting on \k, s) 
the Wigner transformation is W(L(p),k) = [L^ 1 (L(p)k)]L(p)L(k) = L^ 1 {p)L{p)I = /, which consistently produces 
no rotation in spin space, as Ea. l|2.1|l demands. 



III. EXPLICIT GENERATORS OF THE UNITARY REPRESENTATION 



One-particle states arc defined through the action of creation operator a'(p, s) on the vacuum as |p, s) = a^(p, s)|0). 
It follows from Ea. l|2.2|l that under a Lorentz transformation annihilation and creation operators in quantum field 
theory behave as 

a' f (p, S ) = [7(A)at(p, S )C/t(A) (3.1) 
= Jd 3 P [C/* (P)S )(P', S ')( A )]aVy) (3-2) 

= J^D sls (W(A,p))aHp A ,s'), (3.3) 

with p l A = A 1 ^j/ 1 ; and by taking the adjoint, 

a'(p,s) - U(A)a(p,s)U\A) (3.4) 

d 3 P [C/ (P)S)(p ', s0 (A)]a(p )S ') (3.5) 



^^(lU- 1 (A,p))a(p A , S '). (3.6) 



We may hence deduce from the last equation that 



^(p, s )(p',,<)( A ) = y { -^rD ss , (W- 1 (A,p))«(p / - Pa), (3.7) 

and verify by direct calculation, using D^(W) = [DiW)}- 1 = DiW- 1 ), that the transformation U(A) is indeed 
unitary i.e. 

"d 3 p" C/ (p , s)(p » ! ^ ) (A)C/ ( t p „ s „ )(p , s , ) (A) = (5(p-p')^. (3.8) 

But with the formula of the infinitesimal Wigner rotation of Eq . I|2.7fl at hand, it is possible to proceed even further, 
to obtain the explicit generators. By considering an infinitesimal transformation with U (A) = I + ^lo^ li MJ 11 ' we may 
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express Eq. (|3.4|) as 



[/(A)a(p, S )E/t(A) 
= a(p,s) + ^ flv [M^,a(p,s)} 



= (I + w \-^)(l_z[^+ Wo 



■2pf 



k C V 



UjkP j 



lLU 0i 



tjjkJ j p k 

p° + m 



p° + m 

) g + g| 
z<9p l 2ip° 



a(p,s') 



(3.9) 
(3.10) 

(3.11) 
(3.12) 



a(p,s) 



3 „' m/rM" 



d J p 



(p,s)(p',s') 



(3.13) 



In terms of creation and annihilation operators, the explicit infinite- dimensional Hermitian generators of unitary 
transformations U(A) — exp(^o; Ml/ M M,y ) for the non-compact Lorentz group are 



e ijk M jk 



J d 3 pa^{p,s)^J i + ei jk p j j^j a{p,s'), 



= / d 3 paT(p,s) 



jPk 



d 



p u + m 



P 



idp % 2ip { 



p_ 

,-r.O 



a(p,s')- 



(3.14) 
(3.15) 



In general, we may also introduce the particle species label ni for the creation cv(p, s, rii) and annihilation a(p, s, rij) 
operators; and the expression of the generators will then include summing over all rij. 
It can be verified that the explicit generators 



d 

J = J + p x — — 
idp 

j± , J x p n d p . 
p u + m idp 2ip u 

do satisfy the commutation relations of the Lie algebra of the Lorentz group: 



(3.16) 
(3.17) 

(3.18) 



Likewise M.^ of Eqs.(3.14) and (3.15) o bey the similar commutation relations. The expression of the boost generator 
of Ref.[l4l| (following Ref . |20j ) and Ref.|2g which differs from ours in not having the final term of Ea. l|3.17|l can be 
rendered Hermitian provided the measure d 3 p/p° is adopted instead of Weinberg's. 



IV. BELL STATES AND LORENTZ TRANSFORMATIONS 



In this section we specialize to the case with J = cr/2 for spin i particles. For each massive spin | particle, the Lie 
group of all Wigner rotations is SU (2). A two-particle state, and in particular an entangled Bell pair, should transform 
according to the SU(2) x SU(2) representation. As we shall see, the isomorphism between [SU(2) x SU(2)]/ Z2 and 
50(4) permits us to describe the behaviour of the four basis Bell states of spin space for fixed momenta under Lorentz 
transformations succinctly. Since 50(3) is a subgroup of 50(4), in group theoretical terms if we denote the spin i 
doublet as the 2 of SU(2), then a 2-particle state behaves as 

2(g)2 = 4 = l®3 (4.1) 

where the final step denotes its behaviour under 50(3). In other words, we may express the two-particle state in 
terms of a four- vector (the 4) of 50(4) which transforms as singlet and triplet states under 50(3). Indeed it is known 
the four Bell states are expressible as a singlet and a triplet under ordinary 3-dimensional 50(3) rotations, and as we 
shall show, they undergo 50(4) Wigner rotations among themselves under Lorentz transformations. 
The quadruplet (/1 = 0, 1, 2, 3) of Bell states can be conveniently defined in the following manner: 

|B"(p 1)P2 )) = -^ffV^CPl. s ; "i)« f (P2, n 2 )\0) (4.2) 



5 



with <t° = il 2 and a z (i = 1,2,3) being the Pauli matrices a" 1 . In discussing an entangled pair of identical as well as 
distinguishable particles, the additional species labels, n 12 , can be introduced for generality. Let us proceed to show 
that the states defined above are indeed Bell states. The spin indices s, s' are summed over ±5 (which we shall denote 
as ± for simplicity). Focusing on the spin part of the states, (and ignoring for the moment normalization factors and 
species labels), it can be checked that the quadruplet is simply |27j 

15°) ex |+,l;-,2)-|- l;+,2) 
IB 1 }* i(|+,l ;+ ,2)-|-,l;-,2» 
|£? 2 )cx |+,l;+,2) + |-,l;-,2) 
\B 3 ) oc-i(|+,l;-, 2) + |- l;+,2». 

For ease of comparison we have also simplified the momenta indices p 1 2 to 1,2 respectively. Therefore, apart 
from multiplicative constants, these are, respectively, the familiar "singlet" (\B°)) and "triplet" (IB 1 ' 2,3 )) Bell states 
encountered in non-relativistic quantum information science. The conventional assignment of the four Bell states |28[ 
is related to the present one by |/3 U ) = \B°),\/3 W ) = -i\B x ), \j3 00 ) = \B 2 ), |/3 i) = i\B 3 ). But we would like to 
advocate the convention in Ea. (|4.2J) because it is the particular combination of (<7 M a 2 ) ss ia) (p 1; s; ni)a^(p 2 , s'; ri2)\0) 
which provides us with complete and concise description of the behaviour of Bell states under arbitrary Lorentz 
transformations. The upshot is that under generic A, these Bell states undergo a rotation among themselves, and 
transform as 



|g"( Pl ,P 2 ))' = l/(A)|^( Pl ,p 2 )) = ,/ (Api) ° (Ap2) ° fl„"(A, Pl ,p 2 )|g"(p 1A ,p 2A )), (4.4) 

V P1P2 

for which R v M (A,pi,p 2 ) S 50(4) is explicitly listed in Appendix B. 



A. Bell states, and the isomorphism between [SU(2) x SU(2)]/Z 2 and SO(4) 

We shall briefly recap the group isomorphism to establish the notations, and to explain our line of reasoning. Let 
us label elements of the two distinct SU(2) groups by Wi j2 (A). For our purposes, and following the analysis of Wigner 
rotations performed earlier, we should explicitly use 

Wi (A) = exp(i~ • w {Pi)) 
U 2 (A) = exp(i- • w (p 2 )), 

for which 0w(Pi 2 ) are the Wigner angles for particles with momenta p 1 2 . To set up the group isomorphism, we may 
consider 

X = xji» = ( ix ° + x * Xl ~ ix A = ( \ w \ . (4.5) 

M \^El + 1X2 IXq — Xs J \W —V J 

Given U\. 2 G SU{2), it can be verified that X' = IA\XU^~ X is also of the form 

It follows that we may write X' = x'^a^. Since det X — —x^x^ — det X' = —x'x', this implies for each x^ there is 
an R" E SO (A) such that x' — R"x v . Hence we infer 

X' = U x {x^)U^ x = R^x v a", (4.7) 

yielding the identity 

U^U^ 1 = Rfa u . (4.8) 

Returning to one-particle states, the Lorentz transformation for each party of the entangled pair is (disregarding, for 
the moment, normalization factors which are not relevant to the discussion below) 

|Pi,s,ni) 1 — > Wi(A) s / s |p 1A ,s',m) 
|p 2 ,s,n 2 ) 1 — > W 2 (A) s / s |p 2A ,s',n 2 ). 
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The crucial observation is that the quadruplet (/i = 0, 1, 2, 3) of states defined by (a^a 2 ) ss i |p l5 s, m; p 2 , s'; 71.2) trans- 
forms as 

(a'"a 3 ) ss '\p 1 ,s,ni;p 2 ,s',n2) 1 — * [Ui^a 2 )^^ |p 1A , s, m; p 2A , s', n 2 ) (4.9) 

= \U 1 a^U 2 1 o- 2 ] ss ,\p 1A ,s,n 1 \p 2A ,s' ,n 2 ) (4.10) 
= R/^a^^lPiA^^ijpaA.s'^s). (4.11) 

In arriving at the last result we have used a 2 U 2 = U 2 x a 2 which follows from a 2 a 1 a 2 = ~{a l ) T ; as well as the identity 
of Ea. (|4.8|) in the final step. As a consequence, Ea. (|4.4|) is therefore valid. Furthermore Ea. l|4.8|l also yields the 
relation 

RSIA) = ^ ra Tr[Wi(A)a^ 2 - 1 (A)a tt ]. (4.12) 

In passing we mention two special cases: 
For Lorentz transformations which are pure rotations, (X = Y = 9 as 4>(p 1 2 ) = 0); the explicit form of the 5*0(4) 
matrix in Appendix B then yields 



I 



K cos<9% + e ijk (sm 8)6 k + (1 - cos 0)0*0^ 



which means that, as expected, the three Bell states |i3 1 ' 2 ' 3 (p 1 , p 2 )) form an SO(3) rotation triplet while |i?°(p 1 ,p 2 )) 
is a singlet. For an equal-mass entangled pair in the center-of-momentum(COM) frame (p x + p 2 =0) and pure 
boost in the perpendicular direction (t ■ p 1 2 = 0), Eqs. (A9)-(A13) of Appendix A reveal that the resultant Wigner 
rotations are related by 4>{t, p x ) = — 4>(i~, p 2 ) = <p. Thus with (X = Y = <fi) the result is 



(sin 



- (sin <j>) (j) 1 5ij + (cos <f> - 1)4? ft , 



In general, rotational singlet and triplet states do not belong to invariant subspaces when boost transformations are 
also included. Special cases of the behaviour of Bell states for massive particles under Lorentz transformations have 
also been calculated in Refs.0,E3,0]- As we have shown the precise and complete behaviour under arbitrary Lorentz 
transformations can be succinctly stated, as in Ea. l|4.4|l . 

B. Two- particle states as superpositions of Bell states 

In the previous sections we discussed Bell states with infinitely sharp momenta, but it is possible to generalize the 
discussion to generic superpositions of Bell states 

I*) = J d 3 Pl J d 3 p 2 C AI ( Pl ,n 1 ;p 2 ,n 2 )|B' i (p 1 ,n 1 ;p 2 ,n 2 )). (4.13) 

Moreover, it is actually possible to think of any two-particle state, entangled or otherwise, of the form 

\^) = J d 3 p 1 J d 3 p 2 /(p 1 ,si,n 1 ;p 2 ,s 2 ,n 2 )a t (p 1 ,s 1 ,n 1 )a t (p 2 ,s 2 ,n 2 )|0) (4.14) 

in terms of Bell states. Clearly a relation between the coefficients given by 

/(p 1 ,si,ni;p 2 ,s 2 ,n 2 ) = C fl (p 1 ,n 1 ;p 2 ,n 2 )(a tl a 2 ) SlS2 (4.15) 
works. The relation is invertible as 

C M (Pi,ni;p 2 ,n 2 ) = ]^l^f{Pi,s 1 ,n 1 \p 2l s 2l n 2 )((j 2 a u ) S2S1 (4.16) 

Thus C fl (p 1 , m; P2, fi2) can be written down given f(p 1 , si, m; p 2 , s 2 , n 2 ), and vice versa. Note however that in 
quantum field theory all states transform unitarily (1^)' = U(A)\ty)) under Lorentz transformations, no matter how 
complicated the superposition is. The coefficients /(Pi, si, ni] p 2 , s 2 , n 2 ) and C M (p l5 m; p 2 , n 2 ) are not operator- 
valued, and commute with U(A). 
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V. REDUCED DENSITY MATRICES, IDENTICAL PARTICLES, AND LORENTZ-INVARIANCE OF 

VON NEUMANN ENTROPY 

We next introduce reduced density matrices and their properties. While it is possible to generalize, we shall choose 
to concentrate on systems with identical spin | massive fermions e.g. electrons. 

Given an N-particle system of identical particles with density matrix p (which is not restricted to a pure state 
density matrix, but may also correspond to a mixed configuration Trp 2 7^ Trp), the m-particle (m < N) reduced 
density matrices can be defined as 

Pm = ■^^l*i«2---«m)Tr{a il a i2 ...a im pa] m a] m _ l ...aJj(jij2-- Jm|- (5.1) 

This is equivalent to Yang's definition |2l|]. Note that p m is an m-particle operator; and we have simplified all the 
quantum numbers of the creation operator a\ to the label ik- It can then be worked out that Ea. (|5.1|) implies 

(il ■■■i m \Pm \jl---jm) {i\ ■ ■ - hnkm+1 ■ ■ - kn \ Pn \jl ■■■jm^m+l-^n) 



m!(Trp m ) n!(Trp„) 



V to < n , 1 < n < N, (5.2) 



with -£L = ±- . (5.3) 
Trp N Trp 

Thus these reduced density matrices are defined by the partial traces of higher particle number density operators. 

It is worth emphasizing in quantum field theory Lorentz transformations are implemented unitarily on physical 
states. Under any Lorentz transformation A, all creation and annihilation operators transform as U(A)a^ iW (A) and 
U(A)aiW (A). It follows that (we also assume a Lorentz invariant vacuum U(A)\0) = |0)) all states obtained through 
the action of creation and annihilation operators on the vacuum must transform unitarily as \$>}' = C/(A)|\1/) and 
(*S>\' = (\l/|t/T(A). As a consequence, under Lorentz transformations all reduced density matrices defined above also 
transform unitarily as p' = U(A)pU'(A). 

The von Neumann entropy of a density matrix is 

N 

S = -Tr(plnp) = - An In A„. (5.4) 

n=l 

with A„ being the eigenvalues of p (we assume normalization of Trp = 1 has been carried out in the definition of the 
von Neumann entropy). Since p is Hermitian it can be diagonalized, and we may write p = V[diag{X\ 1 A./v)]V^. It 
follows that the eigenvalues are invariant under unitary transformations (and, in particular, Lorentz transformations) 
since p' = U(A)pW (A) = U(A)V[diag(Xi, Xn)][U(A)V]^ obviously has the same eigenvalues as p. Thus the von 
Neumann entropy is Lorentz invariant. Moreover the von Neumann entropy of all reduced density matrices defined 
by S m = — Tr(p m lnp m ) with (m = 1, ...,N) are also invariant for the same reasons. A physical system defined by p 
can thus be parametrized by a set of Lorentz invariant measures {Si, S m , Spj = S}. We shall proceed to show 
that S m can be useful measures of entanglement shortly. 



A. Reduced density matrix, Bell states, and Lorentz-invariant entanglement: a worked example 

Consider a system of two identical fermions, and for the first part of the illustration let us follow the discussion of 
Ref. [25| . A two-particle state of identical fermions may be written as 

I*) = C«ojo}|0) <•../ .A'. (5.5) 

Cij is anti-symmetric and can be set into block diagonal form [2 1) via 

^ for Af odd 

UC U T = (_° C j) 0^ ; (5.6) 

with Aff = (A/"/2) for even Af, and for odd Af, Aff = (Af — l)/2; and U is a unitary matrix. Considering a redefinition 
with a'} = U*j Oj , we may also rewrite 

|*>=2^c ia ^ t +1 |0). (5.7) 

i=i 
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This is the analog of "Schmidt decomposition" for identical fermion systems |23j. The total system has density matrix 
p = and entropy S = 0. Following the prescription for density matrix reduction discussed earlier, the one- 

particle reduced density matrix is then 



(i\ Pl \j) = Tr{a iP a]} = 4(CC%. 
In terms of (UCU T ) combination of Eq.JSJ, we note that Pl = AU^[{UCU T ){UCU T )^]{U)\ thus 

for M odd 



The normalization condition is 



4$>i 2 = i 



Trpi-8^ 



= 2; 



i=i 



and the von Neumann entropy of the reduced density matrix is therefore 

AT, 



Si 



-Tr(^ln^ 
Trpi Trpi 



) = -4]TM 2 ln(2M 2 ); 

8=1 



(5.8) 



(5.9) 



(5.10) 



(5.11) 



which is bounded by In 2 < Si < ln(2A//). The upper limit is obtained by maximizing Si subject to the constraint 
(|5.1UI) : while the lower bound occurs when there is only a single non- vanishing anti-symmetric block in Ea. H5.6fl 23,l25|. 
It may appear disconcerting that this lower bound for fermions is In 2, rather than zero, whereas the analogous result 
for a system of two bosons is < Si < lnJ\f\2Z% . However, the lowest entropy is associated with what we would call 
an "unentangled" system; while any of the "entangled" Bell state discussed earlier will be shown to have Si = 2 In 2 
instead. To see this, let us use, for instance, an unentangled two-fermion state, 



= o+( Pl ,+)a+(p 2 ,-)|0) 
= C l3 a\a)\Q)- 



yielding the non- vanishing part of the C-matrix as C 



(Pir 
(P 2 >- 



(Pi,+) 




P 2 > 



obtain the normalized one-particle reduced density matrix as 



Pi 



* ? 
h 



(5.12) 

(5.13) 
(5.14) 



Following earlier derivations, we 



(5.15) 



giving Si = -2(i) In | = In 2, which is the minimum von Neumann entropy of the reduced 1-particle density matrix 
for two identical fermions. We chose to compute C explicitly to illustrate the consistency of the approach, but it is 
not necessary to go through this step. Given p, computation of pi can also be done directly through Ea. (|5.2l) . 
Consider instead the Bell state 



I*) 



|S°(Pi,p 2 )) 
1 



V2 



a t (p 1 ,+)a t (p 2 ,-) - a t (p 1 ,-)a t (p 2 ,+ 



= £Vt„t 



"!aj-|0>; 



(5.16) 
(5.17) 

(5.18) 



now yielding the non- vanishing part of the C-matrix as 
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(Pn 



(P2>-) (P2>+) (Pi, 



c = 



(Pi. 


+) 


/ 





1 








\ 


(P2> 


-) 




i 













(P2> 


+) 













1 

2\/2 




(Pi. 


-) 


V 








1 

2V2 





/ 



Pi 



, for which the normalized reduced density matrix is 



t\ ' 

I 
\ 

Vo y 



(5.19) 



The corresponding entropy is now Si = — 4(| In |) = 2 In 2 instead. This is a von Neumann entanglement entropy 
(of the reduced density matrix) which is In 2 units higher than the lowest value thus justifying that the Bell state 
|£>°(p 1 ,p 2 )) is an entangled state. Similarly all the Bell states \B fJ, (p 1 ,p 2 )) discussed previously have von Neumann 
entanglement entropy of value 2 In 2. Moreover, as explained, for any physical system the von Neumann entropy of 
any of the reduced density matrices in quantum field theory will be Lorentz invariant. 



B. Generalized zeta functions and von Neumann Entropy 



Zcta function regularizations have been employed in quantum field theory as gauge and Lorentz-invariant methods 
for taming and isolating divergences |29l l30l | . Here we briefly touch upon its relation to the von Neumann entropy. 
The generalized zeta function of an operator can be defined to be 



n s 



(5.20) 



where o„ are eigenvalues of the operator O. There is an interesting relation between generalized zeta function and 
von Neumann entropy. Recall that a density matrix p (or, for this matter, p m , a reduced density matrix) has von 
Neumann entropy 



S = -Tr{p In p) = - } A„ In A„ 



(5.21) 



with A„ being the eigenvalues of p. Note that in the sum, zero eigenvalues do not pose a problem (by L'Hospital 
rule A„lnA n has no contribution for zero eigenvalues). However in quantum field theory, without regularization, 
the entropy can suffer from divergences. For instance, p = e~@ H /Tr(exp"^ ff ) leads to the thermodynamic relation 
S — (3(H), and the expectation value is in general divergent in quantum field theory. 



If we adopt O = p, then ( p (s) = ^2 n jr = ^2 n exp(— s In A„). We note that 



ds 



= — lira \ cxp(— s In A n ) In A n 

s — — 1 ^ — ' 

n 

= A n lnA». 



(5.22) 
(5.23) 



Thus it is possible to define S(s) = -§^-{s), and analytically continue from values for which it is defined to s = — 1. 
However, depending on the form of the operator p, it may still be that s = — 1 is a pole of Cp( s )- 

An alternative then is to note the sum for S in Ea. (|5.21|l does not include zero eigenvalues. Therefore we may 
substitute p with its sub-matrix p' which contains no zero eigenvalues. Its inverse, p' , exists; and we may construct 
its generalized zeta function as 

C[p<]-(s)=I>«) S - ( 5 ' 24 ) 

n 

Since p is a density matrix, its eigenvalues satisfy < A n < 1 Vn; so for positive s, the sum is bounded above by the 
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dimension of p, and therefore converges for any finite-dimensional p. Furthermore 

_ d( p ,-i(s) 



S(s) 



s=l 



ds 



s=l 



lim y exp(s In A„) In A r , 

s— 1 — d 

n 

y^A w ln A„, 



(5.25) 
(5.26) 



and for infinite-dimensional p we may also define the von Neumann entropy S by analytic continuation of the zeta 
function Cp'- 1 ( s ) an d S( s ) to s = 1. A further generalization is 



ds 



(5.27) 



^ =— 1 /a 



with p' substituting for p when a < 0. 



VI. FURTHER COMMENTS AND CONCLUSIONS 



Our work links the Wigner rotations of spins to the behaviour of Bell States under arbitrary Lorentz transfor- 
mations. We discussed reduced density matrices for identical particle systems, established the Lorentz invariance 
of their von Neumann entropies, and suggested an invariant regularization through generalized zeta function. In 
addition, we worked out the explicit expressions of the infinite-dimensional Hermitian generators in the momentum 
representation. We hope that the results presented here will help to place Relativistic Quantum Information Science 
on the surer foundations of quantum field theory which is fully compatible with Lorentz symmetry and causality. It 
is worth emphasizing that we considered the full reduced density matrix with all the degrees of freedom, including 
the momentum with the reduction prescription of Section V when we go from N particles to m < N particles. This 
is different from the reduced density matrix used, for instance in Ref.Q], in which "reduction" to the spin degree 
of freedom is applied even to a single particle. The von Neumann entropy is invariant in our case because Lorentz 
symmetry is unitarily implemented in quantum field theory; each creation operator transforms unitarily under the 
Lorentz group, hence the prescribed reduction of density matrices in Section V will result in Lorentz-invariant von 
Neumann entropy for whatever resultant reduced density matrix we have. In contradistinction, Lorentz symmetry is 
not implemented unitarily in non-relativistic treatments of "wavefunctions" . 

It is known through earlier efforts by others that for a system of two particles in a total pure state, the von 
Neumann entropy of the reduced density matrix is a good measure of the entanglement. We demonstrated that such 
a characterization is in fact Lorentz invariant. Although a system of n particles can be similarly parametrized by the 
Lorentz-invariant von Neumann entropy each of the reduced density matrices, the characterization of entanglement 
in terms of these n numbers, and especially in the case when the total system is a mixed configuration rather than a 
pure state, is still not completely clear and merits further studies. 
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APPENDIX A: COMPUTATION OF THE COMPLETE WIGNER ANGLE 



We shall first calculate the infinitesimal Wigner angle by Taylor expansion as 

u af3 dW 



W(A,p) = W(A,p) 



LJ = 

(A=l) 

uj al 3 r <£L -1 (Ap) 



(A=l) 



a/3 



2 du) a p 



j=0 



dw. 



a(3 



j=0 



^[ L -i(A P )^M 

2 diOa/3 



J a p T -i 



L-\p)M^L(p). 



(Al) 

(A2) 
(A3) 



A=l 



to first order in to; with A(w) = I + ^lo^M^" and A M „ = p^ 1 + Kw^M"' 3 )^ = *t + Straightforward and 
careful calculations using the explicit matrix elements, 



{J )ab i^iab: 



(A4) 
(A5) 



of the generators K l = M 0z and J 1 — \tijkM^ k , and the expression of the Lorentz boost of Eqs. (2.4-2.6) yield the 
complete infinitesimal Wigner rotation as 



W(A,p) = I 



1 



p° + m 



{piLOjO - PjLOio) 



(A6) 



with p° = ^|p| 2 + : 



Ea. i|2.3[l implies Wigner transformations satisfy the multiplication rule 

W(A a Ai,p) = W(A 2 , Kip) ■ W{Ax,p). 



(A7) 



To compute the Wigner rotations for general Lorentz transformations with 6 independent parameters, we may consider 
the decomposition A = R(tp) ■ L(oc). The multiplication rule leads to 



W{A,p) = W(R{ip),L(a)p) ■ W(L(a),p). 



(A8) 



Regardless of the momenta, Wigner angles are degenerate with ordinary rotation angles when the boost parameters 
are zero; thus the first factor is just exp(iip ■ J), while the second factor is the Wigner rotation for an arbitrary 
pure boost. Exploiting the homomorphism between SL(2, C) and 5*0(3, 1), the Wigner angle of this remaining factor 
has been successfully computed by Halpern[3lJ by direct multiplication of 2 x 2 SL(2, C) matrices representing the 
Lorentz transformations in Eq. (2.4-2.6). We may express the result of Halpern as 



with 



W(L(r),p) = exp(z0(r) • J), 

[coshr + coshx + sinhr sinhx(T • p) + (coshr — l)(coshx — 1)(t • p) 2 ] 

[1 + coshr cosh \ + sinhr sinh x(t • p)] 
[to coshr + p° + sinh r(r • p) + (coshr — l)(p° — m)(-f ■ p) 2 ] 

[to + p° cosh r + sinh r(r • p)] 
[sinh r sinh \ + (coshr — l)(coshx — 1)(t • p)] 

[1 + coshr cosh \ + sinhr smhx(T • p)] 
[|p| sinhr + (p° — m)(coshr — 1)(t • p 



(t x p) 



[to + p° cosh r + p sinh r(r • p)] 



(T X p); 



(A9) 

(A10) 
(All) 
(A12) 
(A13) 



and the rapidity \ is related to p by sinh \ = ^ , cosh \ = 

0(r) — > 0! 



It is easy to confirm the infinitesimal limit is indeed 

P X T 



p° + TO 



(A14) 
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The complete expression of the Wigner rotation of (A12) is therefore 

W(A,p) = exp(i0w • J) = exp(iip ■ J) ■ cxp(i</>(a) • J); (A15) 
which yields the explicit relations 

cos(^) = (cos~)(cos~) - (sin-) (sin-) 00 • <£), (A16) 
sm(^-)e w = (cos|)(sin-)0+ (sin-)(cos-)0 + (sin-)(sin-)(0 X 0). (A17) 



APPENDIX B: EXPLICIT FORM OF THE ROTATION MATRIX Rf{A, Pl ,p 2 ) 



The rotation matrix among the four Bell states has been shown to be 



Rf = 7T ^aTr[W 1 (A)a^ 2 - 1 (A)a Q ] (Bl) 



1 

2 

i7? ra Tr[exp(zX • |)^exp(<Y • (B2) 



with the definitions 



X = 9 w ( Pl ), (B3) 
Y = -9 w (p 2 ). (B4) 

(B5) 

Straightforward computations yield the explicit matrix elements 

R ° = (cos|)(cos|)-l-f(sin|)(sin|) 

X Y X Y X Y 

Ri° = -R°i = -(cos — )(sin -)Y t - (sin -)(cos -)X, + (sin -)(sin -)e l]k X 3 Y k 

X Y X Y X Y - X Y 

R. 3 = (cos y)(cos —)5ij - (cos — )(sin — )e iim Y TO + (sin y)(cos —)e ljm X m + (sin y)(sin y)[A- Y5 lj - XiYj - XjYi]. 
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